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ABSTRACT 

We report on 2MASS J01542930+0053266, a faint eclipsing system composed of two 
M dwarfs. The variability of this system was originally discovered during a pilot study 
of the 2MASS Calibration Point Source Working Database. Additional photometry 
from the Sloan Digital Sky Survey yields an 8-passband lightcurve, from which we 
derive an orbital period of 2. 6390157±0. 0000016 days. Spectroscopic followup confirms 
our photometric classification of the system, which is likely composed of M0 and Ml 
dwarfs. Radial velocity measurements allow us to derive the masses (Mi = 0.66 ± 
0.03M Q ; M 2 = 0.62±0.03M Q ) and radii (Ri = O.64±O.O8i? ; R 2 = 0.61 ± O.O9i? ) of 
the components, which are consistent with empirical mass-radius relationships for low- 
mass stars in binary systems. We perform Monte Carlo simulations of the lightcurves 
which allow us to uncover complicated degeneracies between the system parameters. 
Both stars show evidence of Ha emission, something not common in early-type M 
dwarfs. This suggests that binarity may influence the magnetic activity properties of 
low-mass stars; activity in the binary may persist long after the dynamos in their 
isolated counterparts have decayed, yielding a new potential foreground of flaring 
activity for next generation variability surveys. 

Key words: binaries: eclipsing — stars: individual (2MASS J01542930+0053266) - 
stars: low-mass, brown dwarfs 



1 INTRODUCTION 

Low-mass dwarfs (0.07M© < M* < 0.7M Q ) comprise ~ 
75% of all stars in the Milky Way, making them the most 
common luminous objects in the Galaxy (Reid et al. 1995). 
A significant amount of theoretical work has been devoted to 
constructing models that describe the physical processes in 
the interiors and atmospheres of these low-mass stars (e.g. 
Burrows et al. 1993; Baraffe et al. 1998; Hauschildt et al. 
1999; Chabrier & Baraffe 2000). Differences in the predic- 
tions of these models can be subtle, and distinguishing be- 



tween them requires precise empirical constraints on fun- 
damental stellar properties (mass, radius, luminosity, and 
effective temperature; Ribas 2006). As low-mass stars are 
faint, small, and possess complex spectra dominated by 
strong molecular bands, measuring these parameters is chal- 
lenging. Double-lined eclipsing binaries with detached, low- 
mass components of similar spectral type offer the best op- 
portunity for accurate and precise measurements of these 
fundamental properties. 

Although binary systems are more common than single 
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stars at high stellar masses (> 1 Mq; eg. Abt 1983; Duquen- 
noy et al. 1991), binaries are not very common in low-mass 
stars (Leinert et al. 1997; Reid & Gizis 1997; Delfosse et al. 
2004). In fact, the combination of low binary fraction and 
the sheer dominance (by number) of low-mass stars sug- 
gests that most of the stars in the Galaxy are actually single 
stars (Lada 2006). Because of the low binary fractions and 
the faintness of M dwarfs, relatively few low-mass main se- 
quence binaries have been studied (e.g. Creevey et al. 2005; 
Lopez-Morales & Ribas 2005; Hebb et al. 2006; Ribas 2003, 
2006; Blake et al. 2007, and references therein). Adding to 
the census of double-lined eclipsing binary systems is critical 
because they provide highly accurate, direct measurements 
of the masses and radii of the components, nearly indepen- 
dent of any model assumptions. While the known number 
of eclipsing low-mass binaries is small, next-generation vari- 
ability and planet hunting surveys should increase the num- 
ber of known systems. 

Measurements of the masses and radii of the individ- 
ual components of these low-mass systems have revealed 
potentially serious inadequacies in stellar evolution models 
(Ribas 2003, 2006), whereas measurements and models agree 
for stars over IMq. These discrepancies apparently extend 
down into the brown dwarf regime as evidenced by the first 
L-dwarf binary (Stassun et al. 2006, 2007). In particular, 
the mass-radius, mass-temperature, and mass-luminosity re- 
lationships predicted by stellar theory are inconsistent with 
a large fraction of observed M dwarf components (e.g. Hebb 
et al. 2006). Identifying and characterizing new low-mass 
eclipsing binaries will provide stronger constraints for the- 
oretical models and help reveal the cause of the current 
discrepancies between observed and predicted properties of 
low-mass stars. 

In this manuscript, we present the discov- 
ery of a new double-lined eclipsing binary system, 
2MASS J01542930+0053266. In §2, we outline the photo- 
metric and spectroscopic observations of this binary system. 
Our determination of physical parameters and spectral 
types are detailed in §3. Discussion of the discrepancies 
between empirical and theoretical mass-radius relations is 
outlined in §4, followed by our conclusions in §5. 



2 OBSERVATIONS 
2.1 2MASS Photometry 

The calibration observations of the Two Micron All Sky 
Survey (2MASS; Skrutskie et al. 2006) are a precursor to 
next-generation survey efforts such as LSST (Tyson 2002). 
2MASS observed the entire sky using three-channel cam- 
eras that simultaneously imaged in the J (1.25 /im), H 
(1.65 /im), and K s (2.17 fim) passbands. To allow preci- 
sion cross-calibration of the data, 35 1° x 8/5 calibration 
fields were defined. An available calibration field was im- 
aged approximately hourly over the duration of the survey. 
The total number of epochs obtained for a given calibra- 
tion field ranges from 600 to 3700. The 2MASS Calibration 
Point Source Working Database (Cal-PSWDB) has recently 
been released as part of the 2MASS Extended Mission (Cutri 
et al. 2006) , and contains a wealth of information on tempo- 
ral variability (Plavchan et al. 2007; Becker & Cutri 2008). 



2MASS J01542930+0053266 1 (hereafter J0154) was 
first identified as a periodic variable in a pilot study of ob- 
jects in 2MASS calibration region 90004, which was imaged 
2977 times between July 1997 and November 2000. We ex- 
tracted lightcurves for 5770 individual objects from the en- 
semble of data by clustering measurements taken at differ- 
ent epochs based on their celestial coordinates. For a given 
lightcurve, we used the 2MASS photometric quality flag 
phjqual to reject poor-quality data points (phjjual = D, E, 
F, or X). In total, 965 of these clipped lightcurves had more 
than 100 epochs. We phased all clipped lightcurves using a 
modified version of the Supersmoother algorithm (Riemann 
1994). For each object, the J, H, and K s passbands were 
phased independently, and the best-fit periods were com- 
pared. The phasing of J0154 yielded periods of 1.31951246 
days in all three passbands (we note that Supersmoother 
converged upon an alias of the final period derived for this 
system). The composite colours of the system (J — H — 0.66, 
H — K = 0.19) suggest it consists of early-type dM stars 
(Covey et al. 2007). 



2.2 SDSS Photometry 

J0154 lies in the Sloan Digital Sky Survey's (SDSS; Gunn 
et al. 1998; York et al. 2000; Stoughton et al. 2002; Pier 
et al. 2003; Gunn et al. 2006) Stripe 82, a ~ 300 sq. deg. 
equatorial region that has been imaged repeatedly over the 
course of the SDSS. Stripe 82 is imaged by the SDSS-II Su- 
pernova Survey every other night from September to Decem- 
ber, 2005-2007 (Frieman et al. 2008). The extensive repeat 
imaging of this region of sky has enabled precise photomet- 
ric and astrometric calibration of this Stripe, yielding the 
Stripe 82 Light-Motion Curve Catalogue (LMCC; Bramich 
et al. 2008). We extracted the light-motion curve for J0154 
from the LMCC, which consists of 32 u,g,r,i and 2-band 
(Fukugita et al. 1996; Smith et al. 2002) measurements of 
the system, including SDSS-I observations as far back as 
Sept, 1998 and up to the end of SDSS-II supernova survey 
observations in Dec, 2005. The catalogue reports a proper 
motion vector for the system of \i a = 0.88 ± 2.25 mas yr _1 
and /j,s = —11.19 ± 2.25 mas yr _1 . 

Figure 1 displays the ensemble lightcurve, folded at 
the best-fit period of 2.6390157 days. The lightcurves 
are ordered from top to bottom and left to right 
K s , H, J, z, i, r, g, u. The J, H , and K s data are binned every 
30 points. We note that SDSS has no data during the sec- 
ondary eclipse, resulting in poorly constrained relative opti- 
cal colours for each star. This system is one of the faintest 
known eclipsing low-mass systems (r = 18.3), meaning sub- 
stantial telescope time is required to measure the radial ve- 
locity curve to high precision. 



1 For completeness, we note that the averaged position of this 
object in the Calibration Point Source Working Database is 
01542929+0053272. Because the calibration source astrometry 
has a slight bias relative to the survey, having been derived from 
different astrometric standards, we use the designation of the ob- 
ject from the 2MASS All-Sky Point Source Catalog. 
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Figure 2. Spectrum of J0154 obtained on the ARC 3.5-m 
telescope with the DIS-III spectrograph (0.84A/pixel) on 2005 
November 22 UT. Note the Ho emission near 6563 A, indicating 
chromospheric activity. 

2.3 Spectroscopy 

2.3.1 Apache Point Observatory 

To confirm our M-dwarf classification, we obtained spec- 
troscopic observations of this system using the ARC 3.5- 
m telescope with the Dual-Imaging Spectrograph (DIS-III) 
at Apache Point Observatory (APO) on the nights of 2005 
November 22, 2005 November 28, and 2005 December 04 
UT. We used the HIGH resolution gratings (0.84A/pixel in 
the red; 0.62A/pixel in the blue) and a l."5 slit centered 
at 6800A (red) and 4600A (blue). The chips were binned 
2x 1 to increase the signal-to-noise and windowed from their 
original size of 2048k x 2048k to reduce the readout time be- 
tween exposures. The approximate wavelength coverage is 
~ 1000A in both the red and blue wavelength regions. 

The spectroscopic reductions were performed using 
standard I RAF 2 reduction procedures. The calibration im- 
ages (bias, flat, arc, flux) used to correct each of the individ- 
ual spectra were applied only to images taken on the same 
night. The bias, and flat images were observed at the begin- 
ning or end of each night. A He—Ne-Ar arc lamp spectrum 
was taken after each exposure on the target star. 

A representative spectrum of J0154 is displayed in Fig- 
ure 2. These initial spectra clearly demonstrate the features 
of an M-dwarf system. The presence of Ha in emission at 
6563 A indicates magnetic activity, although no obvious line 
splitting is detected. The broad molecular bands of TiO 
(~ 7050 A) are also readily apparent. 

2 IRAF is written and supported by the IRAF program- 
ming group at the National Optical Astronomy Observatories 
(NOAO) in Tucson, AZ. NOAO is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc. under 
cooperative agreement with the National Science Foundation. 
http://iraf.noao.edu/ 



2.3.2 Magellan 

Two 600 second exposures of J0154 were obtained with the 
Low Dispersion Survey Spectrograph 3 (LDSS3) on the Mag- 
ellan II/Clay telescope on the night of December 29, 2005. 
The VPH Blue grating (0.682 A/pixel; R = 1900) and a 
2" slit were used. The spectra were reduced and calibrated 
employing standard techniques, which include subtracting 
a combined bias, subtracting the overscan correction, flat- 
fielding with quartz lamp observations, cleaning the image 
of cosmic rays, and extracting a region centered on the tar- 
get star. The dispersion solution was derived from an ob- 
servation of a He-Ne-Ar arc lamp. Flux calibrations were 
derived from observations of spectrophotometric standard 
stars from Oke & Gunn (1983). The spectra were corrected 
for the continuum atmospheric extinction using mean ex- 
tinction curves. Telluric lines were removed using a proce- 
dure similar to that of Wade & Home (1988) and Bessell 
(1999). We use the Magellan spectra to estimate the spec- 
tral types of the components in Section 3.2.1. 

2.3.3 Keck 

We were unable to resolve line splitting in either the APO 
or Magellan spectra, and therefore made use of the HIRES 
spectrograph at Keck during observing runs on 2006 Octo- 
ber 13, 2006 December 12 and 2007 January 6. Over the 3 
nights we obtained five spectra at R~50,000 with exposure 
times ranging from 30-40 minutes each. 

The data were reduced using standard IDL routines 
that included order extraction, sky subtraction, and cos- 
mic ray removal. The resulting S/N per pixel (0.06 A/pixel) 
ranged from 6-9. The detection of Ha emission lines in both 
stars due to magnetic activity allowed measurement of their 
radial velocities (we were unable to use cross correlation 
techniques on the lower resolution APO and Magellan spec- 
tra). An illustrative example of the Keck data is shown in 
Figure 3. The Ha emission lines were fit to a double Gaus- 
sian profile using Levenberg-Marquardt least-squares min- 
imization. Epochs and radial velocities for each of the 5 ob- 
servations can be found in Table 1. The radial velocity curve 
derived from these data is shown in Figure 4. 



3 ANALYSIS 

3.1 Modeling The System 

We corrected the times of all observations, as well as 
measured radial velocities, to the solar system barycenter 
(TDB). We analyzed the lightcurve of J0154 using the code 
of Mandel & Agol (2002) . All of the orbital elements of the 
binary were allowed to vary, and we allowed the masses 
and radial velocity of the center of mass of the system to 
vary simultaneously while fitting the parameters describing 
the lightcurve. As both stars are well within their Roche 
radii and the centripetal acceleration at the equator is three 
orders of magnitude smaller than the surface gravity, we 
treated each star as a sphere (and hence ignore gravity dark- 
ening). Also, the stars are sufficiently separated that we can 
ignore reflected light. The one uncertainty in our models 
is the limb-darkening, which should be modest in the in- 
frared (Claret et al. 1995). We find that the assumed limb- 
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Figure 3. Spectrum of J0154 obtained with HIRES (0.06 
A/pixel) on the Keck 10-m telescope on Oct 13, 2006. Note the 
line splitting around Ha at 6563 A, yielding emission peaks cen- 
tered near 6561.5 A and 6564.9 A. 
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Figure 4. Radial velocity curve derived from the Keck data. The 
system velocity is 19.1 ± 1.3 km s — 1 . Velocities are relative to the 
solar system barycenter. 

darkening affects our results very little, so we fix the linear 
limb-darkening coefficients for each star, if 1,2, at the val- 
ues computed by Claret (2000, 2004) for model atmospheres 
with T e ff = 3800 and 3600 for the primary and secondary, 
respectively. We assume log(g) = 4.5, [M/H] = 0, and a 
microturbulent velocity of 2 km s _1 for both, based on the 
PHOENIX atmosphere models of Hauschildt et al. (1999). 
The limb-darkening parameters are given in Table 2. 

Our model contains 25 free parameters, starting with 
the five orbital parameters 3 : the eccentricity and longitude 
of pericenter combine to two parameters, ecosw and esino); 
the inclination, i; the time of primary eclipse To (which can 
be translated into a time of periapse); and the period, P 
(which can be translated into a semi-major axis from the 
total mass of the system) . Four parameters describe the bulk 
stellar properties: 7?i, R2, Mi, and M2. The radial velocity 
of the center of mass of the system is 7 (in km s _1 ). The 

3 The longitude of the ascending node is unconstrained since the 
reference plane is the sky plane and nothing in our data constrains 
the position angle of the binary on the sky. 



fluxes are described by 15 parameters (we hold the it-band 
flux of the second star fixed at zero as the best-fit value 
is negative). For the model fitting, we transformed to the 
set of parameters suggested by Tamuz et al. (2006) which 
have weaker correlations between the transformed parame- 
ters. We found the initial best-fit model using Levenberg- 
Marquardt least-squares non-linear optimization giving a 
best-fit model with x 2 = 11304.4 for 9168 de grees of free- 
dom. We found that the scatter of the data outside of eclipse 
had a gaussian shape, but with a larger scatter than the er- 
rors would warrant. It is possible that this discrepancy is due 
to variability in the stellar fluxes as the data were gathered 
over several years, or that the error bars are simply under- 
estimated, so we added a systematic error in quadrature 
(0.04, 0.03, 0.03, 0.01, 0.02, 0.02, 0.01, 0.26 magnitudes in the 
K a , H, J, z, i, r, g, it-bands, respectively) such that the re- 
duced % 2 °f l ne data outside of eclipse in each band is equal 
to unity, and then refit the entire data set. The resulting \ 2 
of 9253.2 has a formal probability P( X 2 > 9253.2) =26% 
for 9168 degrees of freedom. The best fit parameters for the 
brightness of the system are given in Table 2 and the orbital 
and physical parameters in Table 3. 

As the number of free parameters is large and the pa- 
rameters can be strongly correlated, we ran simulated data 
sets, adding gaussian noise to the best-fit lightcurve and ra- 
dial velocity model values. We ran 10 5 simulations, and for 
each simulation re-ran the fitting routine to derive the best- 
fit parameters from the simulated data. We used the 10 5 
sets of derived parameters simulations to compute the l-cr 
errors on the best-fit model parameters, sorting each param- 
eter and choosing the 68.3% confidence intervals, which are 
given in Tables 2 and 3. We computed a 90% confidence 
upper limits on the it band flux of Star 2 by increasing its 
value from zero, while minimizing over the other system pa- 
rameters, until the change in chi-square was A\ 2 = 2.71. 
Our analysis yields asymmetric error bars for all system pa- 
rameters. We quote a single error bar when the negative and 
positive uncertainties are the same to within 10%. 

The best-fit lightcurve is shown in Figure 1. The total 
flux of the two stars is very well constrained due to the large 
number of data points of high photometric quality. However, 
the individual fluxes are more poorly constrained due to the 
small number of points taken during the eclipses, so that 
uncertainties in the radius, limb-darkening, and inclination 
lead to uncertainties in what fraction of each total stellar 
flux is obscured during primary and secondary eclipse. Thus 
the derived fluxes of each star are strongly anti-correlated, 
which is why the errors on the individual fluxes are much 
larger than the error on the total flux (Table 2). The fit 
indicates that, in the infrared passbands, the eclipses are 
of very similar depths (0.35 vs 0.31 magnitudes in the Ks— 
band), which is the reason that the original Supersmoother 
fit converged on an alias of half the period. 

The allowed mass-radius parameter space for each star 
is shown in Figure 5. Each panel shows the probability distri- 
bution of the mass and radius of each star derived from the 
simulated data sets in units of solar radius and mass. The 
contours are l-,2-, and 3-sigma confidence regions (68.3%, 
95.4%, and 99.73% of the 10 s parameter sets). We compare 
the derived values to the masses and radii found for other 
low-mass stars in Section 4.1. 

From our results, we can derive several auxiliary pa- 
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Figure 5. Mass— radius parameter space allowed by our data on J0154. The left figure is for the heavier (primary) object, while the right 
figure is for the secondary. The shaded area correspond to the probability distribution derived from the simulated data sets, while the 
contours represent the l-,2-, and 3-sigma confidence regions. 



rameters for the stars which have different error bars than 
if they are computed from the parameters in Table 3 due to 
covariance between model parameters. We also list in Ta- 
ble 3 (below the horizontal line) : the mass-radius ratio for 
each star, Mip/Rip', the ratio of the radii of the two stars, 
R1/R2; the sum of the stellar radii, Ri + R2; the ratios of 
the semi-major axis to the stellar radii, a/Ri^', the veloc- 
ity semi-amplitudes, K\,% and the total amplitude K; the 
surface gravities, log(gi^)\ the stellar densities p\,2\ and the 
total mass of the system, M1 + M2, which can be used to de- 
rive a semi-major axis of the system, a. At mid-eclipse, the 
projected separation of the stellar centers on the sky is b. 
The fractional error on (Ri +Ri)/a is much smaller (~ 2%) 
than on Ri or R2 individually (> 10%) due to strong cor- 
relations between Ri/a, ife/a and i, as discussed by Tamuz 
et al. (2006). Since our derived fluxes of the stars cover the 
peak in their spectral energy distributions, we have derived 
the bolometric flux ratio of the stars by smoothly interpo- 
lating between the fluxes in the different bands, and then 
taking the ratio of the two stars. Given the relative sizes 
and fluxes of the stars, we derive the ratio of the effective 
temperatures, T2/T1 (we cannot derive the absolute fluxes 
or effective temperatures from the lightcurve and RV data 
as we do not know the distance to this system to high pre- 
cision). 



3.2 Spectral Types and Temperature Estimates 

While the masses, radii and fluxes for each star are well- 
determined by our modeling procedure, there are other 
quantities that can be measured from our observations. The 
spectral types of each component, total space velocity of 
the system and effective temperature estimates can be de- 
termined from our assembled dataset. 



3.2.1 Binary Spectral Template Matching 

In order to estimate the spectral types of this composite 
system, we constructed a grid of binary spectral templates. 
The spectra employed in synthesizing the binary templates 
were drawn from the low-mass stellar templates of Bochan- 
ski et al. (2007) and the K5 and K7 templates used in the 
HAMMER spectral analysis software package (Covey et al. 
2007). Each template spectrum was scaled by its bolomet- 
ric luminosity (Reid & Hawley 2005) and coadded with all 
other templates. Relative velocity shifts were introduced for 
each spectral type pair, ranging from -200 km s -1 to 200 
km s _1 in steps of 20 km s _1 . The final binary spectra grid 
consisted of 1638 templates. The templates were normalized 
to the Magellan observations and residuals were computed 
from 4500 to 7000 A. The best fit binary pair for the Mag- 
ellan spectra is an M0 primary and Ml secondary, with an 
uncertainty of ±1 subtype for each component. We did not 
use the velocity information as a constraint, but rather in- 
cluded the relative shifts for completeness. We adopt these 
subtypes for each component. The best fit composite spec- 
trum is shown in Figure 6, along with the Magellan data. 



3.2.2 Spectral Types From Optical-IR Colours 

Using the colour-spectral type relationships for r — i, i — z, 
and i — J derived by West et al. (2005) and Bochanski et al. 
(2007), we estimate a spectral type of M0 (± 1 subtype) for 
the primary and M3 for the secondary. The intrinsic spread 
in colour at a given spectral type and our measured errors 
permit a range of spectral types between M0 and M4 for the 
secondary. Because of this large uncertainty, we adopt the 
secondary subtype derived from the spectral template anal- 
ysis, and note that these colour-based results are consistent 
with the spectroscopic results derived in Section 3.2.1. We 
estimate a distance of 623 ± 50 pc from the i-band photo- 
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Figure 6. Magellan spectra (0.682 A/pixel) of J0154 along with 
the best fit binary spectral template. The best fit binary template 
is composed of an M0 + Ml. 



metric parallax (West et al. 2005) of the primary star. At a 
Galactic latitude of b — —58°, this binary has vertical dis- 
tance below the Disk of ~ 530 pc, consistent with being a 
member of the Galactic thin disk. 

The thin disk membership of J0154 is strengthened by 
examining the system's kinematics. The tangential veloc- 
ity implied by the observed proper motion (11.2 ± 3.2 mas 
year -1 ) and distance to the system is 33.1 ± 12.1 km s _1 . 
Added in quadrature with the system radial velocity from 
the binary fit, we find a space velocity of 38.2 ± 12.2 km s~\ 
again consistent with the thin disk (Bochanski et al. 2007). 



3.2.3 Metallicity and Temperature Estimates 

We measured the composite CaH2, CaH3, and Ti05 molecu- 
lar indices in our APO spectra of this system, and find CaH2 
= 0.64 ±0.06, CaH3 = 0.83 ±0.03, and Ti05 = 0.64 ±0.04. 
Using these values along with the empirical formula from 
Figure 2 of Woolf & Wallerstein (2006) yields an effective 
temperature estimate for the primary of T e s — 3730±100A', 
consistent with the M0 spectral type determined from the 
full spectra. Further, when these measurements are com- 
pared with the samples of Lepine et al. (2003), Woolf & 
Wallerstein (2006), Burgasser & Kirkpatrick (2006), and 
West et al. (2008), they suggest that the composite system is 
of solar or slightly super-solar metallicity, again consistent 
with thin disk membership 4 . 



4 DISCUSSION 

Previous studies (e.g. Ribas 2006, and references therein), 
have demonstrated that current models underpredict the 
radii of low-mass stars at a given mass. The radii derived 
for J0154's components (Ri = 0.64 ± O.O8i? at Mi = 
0.66±0.03M Q , R 2 = 0.61 ±0.09^0 at M 2 = 0.62±0.03M S ) 
lie in between different model predictions, and the errors 
are currently large enough to not yet provide discrimina- 
tion between the models. The source of the uncertainties 

4 It should be noted that there exist dMs in Woolf & Wallerstein 
(2006) with similar CaH and TiO indices to our targets, but with 
[Fe/H] values below -0.3 




Figure 7. Confidence limits on b/(Ri + R2) and R2/R1 from the 
synthetic data sets. The shaded region shows the density of points 
from the best-fit parameters to 10 5 simulated data sets, while the 
contours are 68.3%, 95.4%, and 99.73% confidence regions. 



on the stellar radii is almost entirely due to the large er- 
rors on the photometry given the faintness of this system. 
A severe banana-shaped degeneracy between the impact- 
parameter (or inclination) and the ratio of the stellar radii 
occurs due to the large photometric errors (Figure 7); this 
is why in Table 3 the fractional error on Ri ± R2 (~ 3%) is 
much smaller than the fractional error on R2/R1 (~ 30%). 
Within the 68.3% confidence limit, the deviation of the K- 
band lightcurve from the best-fit lightcurve is only 0.6%, 
which implies that milli-magnitude precision would be re- 
quired to derive the radius ratio to high accuracy. Below, 
we discuss the implications of our system with regards to 
current theoretical and empirical mass-radius relations. 



4.1 The Empirical Mass— Radius Relationship 

In Figure 8, the masses and radii of known low-mass eclips- 
ing binary systems (Lopez-Morales 2007; Lopez-Morales 
et al. 2006; Lopez-Morales & Shaw 2007) are plotted along 
with current empirical (Bayless & Orosz 2006) and theoret- 
ical models (Baraffe et al. 1998; Siess et al. 2000). Bayless 
& Orosz (2006) have derived an empirical mass-radius re- 
lationship for K and M dwarfs from known binaries that 
stretches up to nearly O.8M0. We test this empirical mass- 
radius relation by comparing their expected radii, given our 
mass measurements, to our measured radii. Their analysis 
predicts Ri = O.67±O.O3i? and R 2 = 0.63 ± O.O3i? , while 
we measure Rx = 0.64 ± O.O87? and R 2 = 0.61 ± O.O9i? . 
Our objects are consistent with the ensemble of eclipsing 
binary stars used to derive their relationship, and fill in the 
gap at the high-mass, early-dM end of the relationship. The 
dearth of data in Figure 8 and the recent discovery of many 
of these systems reflects that this is an emerging field, only 
recently enabled by large-scale photometric surveys. 

Theoretical models predict mass-radius relations which 
are a strong functions of both metallicity and age (Baraffe & 
Chabrier 1996). The predictions of the Baraffe et al. (1998) 
evolutionary models for objects of solar metallicity and ages 
of 10 s -10 9 years are consistent with J0154, yet disagree with 
other systems of similar mass. The models tend to system- 
atically underpredict the radii at a given mass, as shown in 
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0.2 0.3 0.4 0.5 0.6 0.7 0.8 

M(M 9 ) 

Figure 8. Distribution of mass vs. radius for known low- 
mass double— lined eclipsing binary systems (open circles; Lopez- 
Morales 2007; Lopez-Morales et al. 2006; Lopez-Morales & Shaw 
2007; Maxted et al. 2004). The solid line shows an empirical rela- 
tionship derived by Bayless &; Orosz (2006); the dotted line rep- 
resents the Z=0.02, Y=0.275, T=l Gyr models of Baraffe et al. 
(1998); the dashed line represents the Z=0.02, Y=0.28, T=l Gyr 
models of Siess et al. (2000). The solid dots correspond to the 
components of J0154. 

Figure 8. While the uncertainties on the masses and radii of 
J0154's components render them an equally good fit to both 
the Baraffe et al. (1998) theoretical models and Bayless & 
Orosz (2006) empirical fit, additional photometric and spec- 
troscopic data should yield more precise measurements of 
these attributes, and better discrimination between models. 

For comparison, the models of Siess et al. (2000) are 
also shown in Figure 8. The large differences between mod- 
els with similar inputs for metallicity and age highlight the 
uncertainty that presently exists in this field. Hopefully this 
situation will be remedied by more high-precision measure- 
ments of fundamental stellar parameters in binary systems, 
along with updated models. 

4.2 Activity 

An interesting aspect of this system is the observed Ha emis- 
sion in both of the components. West et al. (2004) find that 
less than 5% of isolated M0 and Ml stars show activity (Ha 
equivalent width of at least 1A). The activity in early- type 
M dwarfs is also short lived (< 1 Gyr; Hawley et al. 2000; 
West et al. 2008). It would be very unlikely to randomly 
draw an active M0 and Ml star from the field population, 
suggesting that some aspect of the interactions between the 
components is inducing the observed magnetic activity. The 
improbability of the stars being independently active is fur- 
thered when we consider the distance that this pair is from 
the Galactic plane. It is likely that M0 and Ml dwarfs that 
are several hundred pc from the Plane have been dynam- 
ically heated for many Gyrs and have ceased being active 
(West et al. 2008). 

A search through the Palomar/MSU Nearby Star Cata- 
log of Gizis et al. (2002) shows that a large fraction (20/22) 
of double-lined spectroscopic M-dwarf binaries have mag- 
netically active components (Ha equivalent width of at least 
lA). Other empirical studies (e.g. Lopez- Morales 2007) have 



suggested that the magnetic activity and metallicity of a 
star may affect its radius, drawing an explicit connection 
between the enhanced activity and large radii found in M- 
dwarfs in binary systems. Chabrier et al. (2007) suggest that 
enhanced magnetic activity may lead to inefficient thermal 
transport in stellar interiors. This results in objects with 
larger radii and smaller effective temperatures than stars 
where magnetic effects are negligible. Rapid stellar rotation 
may similarly affect the interior convection. In addition, en- 
hanced surface spot coverage (30 — 50%), due to strong mag- 
netic activity, may impact the stellar radius. Either of these 
effects may be responsible for the empirical mass-radius re- 
lationships found in low-mass binary stars; additional data 
are required to constrain these theories. 



5 CONCLUSIONS 

We report the discovery and characterization of the double- 
lined eclipsing binary system 2MASS J01542930+0053266. 
Photometric and spectroscopic evidence suggests that both 
components are M dwarfs, and we adopt classifications of 
M0 and Ml as their subtypes. We resolve splitting of the Ha 
emission line with spectroscopic observations using HIRES 
at Keck, leading to a radial velocity curve and estimates of 
the masses and radii of each star. Simulated data sets cre- 
ated by adding noise to the best fit provides uncertainties 
on and covariances between the system parameters. We em- 
phasize that there exist complicated degeneracies between 
parameters in eclipsing systems that can only be fully ex- 
plored with such detailed analyses. 

Our analysis is consistent with previous studies of 
double-lined eclipsing M-dwarf systems. An empirical study 
by Bayless & Orosz (2006) yields a quadratic mass-radius 
relationship spanning the range of 0.2 — 0.8Mq and spec- 
tral types from late K to late M. This empirical relation 
accurately predicts the radii of J0154's components. 

We observe Ha emission from both components, an un- 
likely scenario given their early-M spectral types and their 
distance from the Galactic plane. If magnetic activity is en- 
hanced in M dwarfs in binary systems, the binary population 
including M-dwarfs components may present an additional 
foreground of stellar flares for next generation time domain 
surveys (Kulkarni & Rau 2006; Becker et al. 2004). 
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Table 1. Radial Velocities 



Date (TDB) 


v\ (km s 


- 1 ) 


t>2 (km s 1 ) 


54021.3758 


-55.49 ± 


2.75 


99.51 ±2.95 


54021.4043 


-57.00 ± 


2.89 


100.90 ± 3.25 


54081.2152 


81.77 ±2.80 


-47.93 ± 3.45 


54106.2401 


-41.12 ± 


3.06 


91.28 ±3.74 


54106.2778 


-46.12 ± 


3.12 


80.38 ±3.38 



Baryccntric radial velocities measured from the Keck data on Oct 
13, 2006, Dec 12, 2006 and Jan 6, 2007. The errors include a sys- 
tematic error of 1.75 km s — 1 which has been added in quadrature 
to each of the measured velocity uncertainties. Dates are in Mod- 
ified Julian Day (MJD) corrected to the solar system dynamic 
barycenter (TDB). 

West A. A., Walkowicz L. M., Hawley S. L., 2005, PASP, 
117, 706 

West A. A., et al. 2008, ApJ, Submitted 

Woolf V. M., Wallerstein G., 2006, PASP, 118, 218 

York D. G., et al., 2000, AJ, 120, 1579 
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Table 2. Apparent brightness, colours, and limb-darkening of J0154 



Filter 



Star 1 



Star 2 



Total 



/2//1 



mi - K s 



m-2 - ^s,2 rntot - K. 



s.tot 



Ul 



«2 



U 

9 
r 



J 
H 



21.95 ±0.05 
19.89+g; 32 
18.59+0^ 



> 22.74" 



17.69 



-0.23 
+0.32 
-0.23 
1 7 09+O.3I 
i( - ZZ -0.23 

16.08+^ 
15.4«±g;« 
15.28+0L 1 



21.26 
19.70 



+ 1.22 

-0.72 

-1.06 

-0.59 

18.66+^2 
17.99+°^ 
16 47+° :42 

10 -*' -0 .32 

15 74+ - 41 
10 - '^-0.30 
15 ,-0+0.40 



21.9477 ±0. 
19.6159 ± 0. 
18.2573 ± 0. 
17.3164 ± 0. 
16.7874 ± 0. 
15.5019 ± 0. 
14.8373 ± 0. 
14.6430 ± 0. 



0461 
0038 
0039 
0035 
0032 
0013 
0015 
0020 



< 0.95" 
0.2833+°-« 
0-3581+^ 

0.4100t°-^ 
0.4943+°j0 

0.69791 0.31 
0.76721°-^ 
0.7971±£«j 



6.67 



+0.247 
-0.323 



4.6V 



■ +0.072 
-0.057 

3.31 ±0.07 
2.41 ±0.07 
1.94 ±0.06 
0.80 ±0.05 
0.18 ±0.05 



> 7.2" 

c 7 o+0.875 

A 1 q+0.670 
4 ' 1S -0.303 
o 1 n+0. 505 
°- lo -0.238 

9 46+ - 298 
Z - 4D -0.167 

0.94 ± 0.06 
0.22 ± 0.06 



7.3047 
4.9729 
3.6143 
2.6734 
2.1444 
0.8589 
0.1943 



± 0.0462 
± 0.0043 
± 0.0043 
± 0.0041 
± 0.0038 
± 0.0024 
±0.0025 



0.713 
0.829 
0.808 
0.696 
0.611 
0.481 
0.453 
0.377 



0.734 
0.814 
0.777 
0.672 
0.589 
0.428 
0.398 
0.329 



a : In the u band the best-fit flux for the second star is zero, so we report 90% limits on the magnitude and colours of this star. 
The apparent brightnesses of each component of J0154, as well as the composite system brightness, as derived from our global fit. We 
also list the adopted limb— darkening parameters U12 for the components (Claret 2000, 2004). The total system brightness is better 
constrained than those of the individual components because of the latter's covariancc with the radii of the stars, as well as the 
inclination of the system. Due to the anticorrelation between the fluxes of the two stars, the flux ratio (column 5) has larger 
uncertainties. The colours have smaller errors than the magnitudes since the fluxes in each band are strongly correlated for each star. 
We quote a single error bar when the positive and negative uncertainties are the same to within 10%. 



Table 3. Binary Fit Parameters 



Parameter 



Value 



e cos (a;) 
i (radians) 
e sin(w) 
T (TDB) 
P (days) 
Ri(Rq) 
R 2 (R @ ) 
Mi(M ) 
M 2 (M Q ) 
7 (km s _1 ) 



0.00142 ± 0.00068 
1 51445+0.00839 

1.01440_p .00109 

-0.006 ±0.012 
52244.82052 ± 0.00058 
2.6390157 ± 0.0000016 

639+ ' 076 

610+ ' 083 

U.D1U_q 1Q2 

0.659 ±0.031 
0.619 ±0.028 
19.09 ± 1.28 



Mi/Ri(M Q /R e ) 

M 2 /R 2 (Mq/R & ) 

i?i + R 2 (Rq) 

R 2 /Ri 

a/R-i 

a/R 2 

K\ (km s- 1 ) 
K 2 (km s- 1 ) 
K (km s- 1 ) 
log(gi [cm s~ 2 ]) 
'09(92 [cm s~ 2 ]) 
pi (g cm~ 3 ) 
p 2 (g cm" 3 ) 
Mi + M 2 (M e ) 
M 2 /Mi 
a (R Q ) 
a (AU) 

b (Rq) 

T 2 /Ti 



1.032Z 
1.016j 
1.248^ 
0.955^ 
13.652^ 
14.299" 1 



-,+0.171 
2 -0.112 
,+0.204 
" -0.124 
j+0.021 
-0.048 
-+0.304 
-0.244 
12.220 
-1.455 
5+2.850 
-1.707 

80.896 ± 1.811 
76.062 ± 4.204 
156.959 ± 5.901 
4 646+ 131 

t.DtD_ Q9g 

4 659+ ' 158 

q cco+2.032 
o.owij ^ 023 

3.849+. 2 ^° 
1.278 ±0.054 
0.940 ±0.034 
8.718 ±0.123 

0.041 ±0.001 

,o 6 +0.005 

U - M ^'.016 



Best fit system and physical parameters of J0154 (above line) and derived parameters (below line). The fitting process is described in 
Section 3.1. Uncertainties in the parameters are derived from the distribution of best- fit values for the 10 5 synthetic light curves. We 
quote a single error bar when the positive and negative uncertainties are the same to within 10%. 



